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Abstract. In recent works, we have proposed a stochastic cellular automaton model 
of traffic flow connecting two exactly solvable stochastic processes, i.e., the Asymmetric 
Simple Exclusion Process and the Zero Range Process, with an additional parameter. 
It is also regarded as an extended version of the Optimal Velocity model, and moreover 
it shows particularly notable properties. In this paper, we report that when taking 
Optimal Velocity function to be a step function, all of the flux-density graph (i.e. the 
fundamental diagram) can be estimated. We first find that the fundamental diagram 
consists of two line segments resembling an inversed-X form, and next identify their 
end-points from a microscopic behaviour of vehicles. It is notable that by using a 
microscopic parameter which indicates a driver's sensitivity to the traffic situation, we 
give an explicit formula for the critical point at which a traffic jam phase arises. We 
also compare these analytical results with those of the Optimal Velocity model, and 
point out the crucial differences between them. 
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1. Introduction 

Traffic dynamics lias naturally attracted much attention from engineers since the volume 
of vehicular traffic outstripped road capacity [H [2|, [3]. Especially for the last decades, 
it has attracted a lot of interest, in addition, from physicists and mathematicians as 
a typical example of non-equilibrium statistical mechanics of self-driven many particle 
systems [1]. Since self-driven particles do not obey Newton's laws of motion, their 
collective phenomena are far from predictable and are strongly dependent on their 
density. Recently, a number of different approaches to the subject have been made 
from various viewpoints such as microscopic and macroscopic, continuous and discrete, 
deterministic and stochastic. 

The Burgers equation pt = 2ppx + Pxx-, is known as an elementary model of 
traffic flow. It formulates the evolution of density distribution of vehicles, i.e., it is 
a macroscopic, continuous model. Recent studies show that the Burgers equation is 
directly connected with other basic models of traffic flow [4j. The Burgers equation is, 
at ffist, transformed into a cellular automaton (CA) model (a microscopic discrete one) 
through ultra- discretization [5]. Note that CA models are microscopic ones because they 
explicitly define each particle's motion. It provides a method to reveal another profile 
of models, and enables us to return from macroscopic to microscopic. From the Burgers 
equation, we thereby obtain the so-called Burgers cellular automata (BCA), which is 
an extension of the Rule-184 CA [4j. This CA has a very simple update rule, i.e., if 
the adjacent site is not occupied, the vehicle move ahead with a given probability, but 
otherwise it does not. Each site contains one vehicle at most, and the rule is generally 
called the hard-core exclusion rule. In the case of BCA, each site can contain more than 
one vehicle. 

Moreover, BCA can be transformed into another basic model, i.e., the Optimal 
Velocity (OV) model, through the discrete Euler-Lagrange (E-L) transformation. The 
discrete E-L transformation is made on fully discrete variables, and then field variables 
change to particle variables |4j. The OV model is a continuous, microscopic model, and 
is expressed as 



where Xi = Xi{t) denotes the position of the i-th vehicle at time t and the function V 
is called the Optimal Velocity (OV) function [6l[7j. The OV function gives the optimal 
velocity of a vehicle in terms of the headway Xj+i — Xi, where the 2-th vehicle follows 
the {i + l)-th in the same lane, and then the OV function, in general, is monotonically 
increasing. In particular, the OV model obtained from BCA has OV function which is 
a step function. This suggests that the OV model with a step function is essential as 
well as elementary. 

Cellular automaton models are efficient and flexible compared to those described 
by differential equations, and they have been used to model complex traffic systems 
such as ramps and crossings [2]. The Nagel-Schreckenberg (N-S) model, a well-known 
CA model, successfully reproduces typical properties of real traffic [8]. What makes it 
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sophisticated is a random braking rule, which is a plausible mechanism to simulate the 
motion of vehicles in a single lane. Nevertheless, the N-S model does not succeed in 
reproducing the so-called metastatic state, i.e., an unstable state which breaks down to 
the lower-flux stable state under some perturbations. 

Extensive study of traffic flow has revealed that the metastable property, appearing 
in the medium density region, is universal in real traffic flows, and accordingly that 
property is essential in modelling [9l |10]. Moreover, it is quite distinct among non- 
equilibrium statistical systems [I]. In other words, this metastable property plays a 
critical role in characterizing traffic flow from the viewpoint of dynamics, and it is hence 
required for traffic models to exhibit this property. Thus far, one needs the slow-start 
rule to reproduce a metastable state in existing CA models [lOl [TTl [121 113] . It introduces 
a delay for vehicles to respond to the changing traffic situation, i.e., if a vehicle stops 
due to the hard-core exclusion rule, the slow-start rule forces it to stop again at the 
next time step. In contrast, due to the second-order derivative, the OV model naturally 
includes a similar mechanism to the slow-start rule. The intrinsic parameter a in the OV 
model ([1]) corresponds to the driver's sensitivity to a traffic situation (e.g. the distances 
relatively to the vehicles ahead), and plays an important role in the stability of a traffic 
flow. Note that the reciprocal of a represents a driver's response time, connecting the 
OV model and the Newell model [H]. 

In the next section, we introduce a stochastic CA model following and then 
we show that the model inherits the sophisticated features from the OV model. 

2. Stochastic optimal velocity model 

Most stochastic models incorporate noise, taking into account uncertain effects such 
as different driver characteristics, driver error false operation, and external influences. 
In general, randomness disturbs metastable states and thus stochastic models do not 
show any metastable state in the fundamental diagram. In contrast, as will be seen 
in the latter part of this paper, the probability of the SOV model plays an essential 
role in producing metastable states, and the metastable states emerge by extracting a 
deterministic mechanism. Note that a similar trick was considered in [2T] . 

2.1. General scheme 

First of all, we explain the general framework of our stochastic CA model for one-lane 
traffic. The roadway, being divided into cells, is regarded as a one- dimensional array of 
L sites, and each site contains one vehicle at most. Let M* be a stochastic variable which 
denotes the number of sites through which the i-th vehicle moves at time t, and wl{m) 
be the probability that M ■ = m (m = 0, 1, 2, . . .). Then, we assume a principle of motion 
that the probability wj'^^{m) depends on the probability distribution wj{0),wj{l), . . ., 
and the positions of vehicles x\, x\, . . . , x^j^ at the previous time. The updating procedure 
is as follows: 



Analytical study on the criticality of the Stochastic Optimal Velocity model 4 

• Calculate the next intention wl~^^ {i = 1,2,..., A^) from the present, intention 
w*(0), . . . and positions x\,x\, . . . , x^; 

wl-^'im) = fiwm,wl{l), . . . ; X* , . . . , x^; m) (2) 

• Determine the number of sites M*"*"^ through which a vehicle moves (i.e. the velocity) 
probabilistically according to the intention wj'^^. 

• The new position of each vehicle is 

= xl + min(Ax*, M*+i) (Vz), (3) 

where Ax* = x^^-^ — x\ — 1 denotes the headway. (Headway is defined to be the 
clear space in front of the vehicle, and thus in a CA model we need to subtract 1 
to take account of the site occupied by the vehicle itself.) 

The hard-core exclusion rule is incorporated through the second term of the right hand 
side of (E]). 

We call the probability distribution w\ the intention because it is an intrinsic 
variable of the vehicle and drives themselves. It brings uncertainty of operation into the 
traffic model and has no physical counterpart. 



2.2. The SOV model 

In what follows, we assume w\{m) = for m > 2. It is notable that Yl,^=Qw\{m) = 1 
by definition and the expectation value (M*) = J2m=o^'^ii^)y hence, setting 
f • = wl{l), we have w*(0) = 1 — f * and (M*) = f*. From ([2]) we have 

/ wt^'il) = vj^' = f^l;xlxl...;r) 

and we therefore express the intention by vj in stead of wj. As long as vehicles move 
separately (i.e. Ax* ^ 0), the positions are updated according to the simple form 

J X* + 1 with probabihty vl~^^ , . 

' \ x\ with probability 1 - f -^^ 

and consequently we have 

(^^^) = (^*) + ^r (6) 

in the sense of expectation value. This equation expresses the fact that the intention 
v\'^^ can be regarded as the average velocity at time t. 
Let us take an evolution equation 

t;*+^ = (l-a)t;* + ar(Ax*), (7) 

in ([2]), where a (0 < a < 1) is a parameter and the function V takes the value in [0, 1] so 
that v\ should be within [0, 1]. Equation ([7]) consists of two terms, i.e., a term turning 
over the intention v\ into the next, and an effect of the situation (the headway Ax-). 
The intrinsic parameter a indicates the sensitivity of vehicles to the traffic situation, 
and the larger a is, the less time a vehicle takes to change its intention. 
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A discrete version of the OV model is expressed as 

Xi{t + At) - Xi{t) = Vi{t)At, 

Viit + At) - Viit) = a\v{Axi{t)) - Vi{t)]At, 



(8) 
(9) 



where Axi{t) = Xi+i{t) — Xi{t), and At is a time interval Due to the formal 
correspondence between ((Tj) and ([9]), we call a stochastic CA model defined by (171) 
the Stochastic Optimal Velocity (SOV) model, hereafter. 

As we noted in the preceding work [15], the SOV model reduces to two exactly 
solvable stochastic models, the Asymmetric Simple Exclusion Process (ASEP) [161 [13 
and the Zero Range Process (ZRP) [iHl [19], when the parameter a (0 < a < 1) takes 
the values of and 1 respectively. The SOV model has a fundamental diagram similar 
to that of ZRP as a approaches 1 (the ZRP limit). However, when a takes a small value 
(the ASEP limit), the SOV model shows quite different properties from those of ASEP. 
Since we already discussed this point in detail |[15], we do not devote any space to the 
limiting case a — > 0. It should be noted here that a model including ZRP and ASEP 
was proposed in |20J. However, this model does not satisfy the essential standard that 
a traffic model should reproduce the metastable state observed universally in empirical 
traffic data. 

3. The SOV model with a step OV function 

In this section, we take a step function 



as the OV function. Then, we impose a periodic boundary condition and adopt the 
parallel updating as usual. 

3.1. Comparison with the OV model 

In preceding works [221 [23], the original OV model with the step OV function ( fTOl) 
has been studied in detail. In an ideal case, they gave an exact solution for cluster 
formation (i.e. the traveling cluster solution). They use figures to illustrate the 
behaviour of vehicles forming a cluster, and find a closed loop or hysteresis loop in the 
velocity-headway diagram (phase space). Moreover, from the microscopic viewpoint, 
they evaluate the velocity of a jam transmitting backward, and two specific headways; 
the distance Axj with which vehicles stop in a jam, and Axp with which vehicles move 
in a free-fiow region. (See §3.4[ ) 

Figure [H shows the phase space of our model. Note that the vertical axis of Figure [H 
means the expectation value of velocity. We find that there appears a hysteresis loop 
corresponding to Figure 3 in [22] . Since our model incorporates the hard-core exclusion 
rule as well as a CA model with maximum allowed velocity 1, it takes a different shape 
from the parallelogram of Figure 3 in [22]. Note that, in contrast with the original OV 





(0 < X < rf) 
{x > d) 



(10) 
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Figure 1. The headway-velocity diagram (phase space) of the SOV model (a = 0.8) 
with a step function (d = 2). We observe a hysteresis loop (thick gray line) around 
the discontinuous point of the step OV function (thin black line), where the phase of 
a vehicle goes round in the direction of arrows. In comparison with Fig. 3 of 22J, two 
corners of the parallelogram are folded by an effect of the hard-core exclusion rule. 



Position 




Figure 2. The spatio-temporal pattern of the SOV model (a = 0.8) with a step 
function {d ~ 2) simulated with the number of site L — 1000 and the number of 
vehicles M — 400, where periodic boundary condition is imposed on the roadway. 
There appear a lot of stable clusters (small jams) propagating backward. 

model, our model contains the element of randomness and hence the vehicles do not 
always move in accordance with the hysteresis loop. In Figure [2], we show the spatio- 
temporal pattern of our model, which should be compared with Figure 4 in [22]. Due to 
the randomness, the clusters do not hold their sizes constant and the number of them 
fluctuates with time. 

These two figures suggest that the SOV model inherits, from the original OV 
model, the mechanism of vehicles clustering and then separating. However, from the 
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viewpoint of many-particle systems, the SOV model shows apparently different collective 
phenomena from that of the OV model. In the following part of this section, we estimate 
the fundamental diagram of the SOV model with a step function as well as comparing 
the estimated specific headways and velocity of our model with those given in [22]. 



3.2. Fundamental diagram 

We denote the density of vehicles to sites by p = N/ L (L is the number of sites, and the 
number of vehicles), which is a macroscopic variable, and a conserved quantity of motion 
under the periodic boundary condition, no entrances or exits. Another macroscopic 
variable flux, Q = pv, is defined using the average velocity in a steady state; 

v:=^U^J-xJ-'), (11) 
1=1 

where time T should be taken large enough for the system to reach a steady state. 

A fundamental diagram, a plot of the flux versus the density, illustrates how traffic 
conditions depend on density. It represents the characteristics of a traffic model, and 
hence traffic models are required to reproduce a fundamental diagram observed in real 
traffic flow. As far as the above-mentioned exactly solvable models (ASEP and ZRP) 
are concerned, we can make an exact calculation of the fundamental diagram [25t [26]. 
We only exhibit an explicit formula of the fundamental diagram of ZRP (i.e. a = 1) 
especially when a step function is adopted as the OV function: 

^ 1 min(p, 1-dp) (0 < P < ^) (12) 

where d is the discontinuous point of the step OV function (ITUi) . 

Since the flux of R184-CA (ASEP with the hopping probability p = 1) is obtained 
in the simple form of Q{p) = min(p, 1 — p), considering that each vehicle occupies d 
sites, we directly obtain the formula (fT2l) . 



3.3. Metastable state 

Figure [3] shows the fundamental diagram of the SOV model with the step OV function 
(fTOj) . We flnd that the diagram consists of two lines corresponding respectively to free- 
flow phase (positive slope) and jam phase (negative slope). It is remarkable that there 
is a region of density where two states (a free-flow state and a jam state) coexist. The 
second-order difference allows the model to show this property, so-called hysteresis, as 
pointed out in Section 1. 

The free-flow line in the fundamental diagram has a slope of 1, i.e., all the vehicles 
are moving deterministically (i.e. vj = 1) without jamming. This kind of state can 
be implemented in the uniform state; equal spacing of vehicles and the initial velocity 
= 1 for all i. As can be seen from ([7]), the intention changes 

vl = l-il-v^){l-ay iAx\>d), 

vl = vfil~ ay (Ax* < d), ^ ' 
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Figure 3. The fundamental diagram of the SOV model with the step OV function 
(|10P plotted at each value of sensitivity parameter a, where the discontinuous point is 
d = 2, the initial value of the intention is vf — 1, and the system size is L = 1000. 
Theoretical curve (gray) has a complete agreement with simulated data (dots) in these 
cases. 



(a) a=0.2 



(b) a=0.8 



1^ 



X 



500 1000 1500 

Time 



1000 1500 2000 

Time 



Figure 4. The flux decreases rapidly as an external perturbation is imposed, and 
finally settles into a lower value, i.e., the flux of jam state. We simulate with the number 
of sites L = 1000, the number of vehicles M = 334, and the sensitivity parameter takes 
the value of (a) a = 0.2 and (b) a = 0.8. It imply that, in contrast with our preceding 
paper, the metastable state has no lifetime. 



and consequently the uniform states constitute a line segment 

Q = P (0 < p < Ph) (14) 

in the fundamental diagram. The maximum- flux density ph, at which the vehicles take 
the minimum value d of equal spacing, is given as follows: 

P. = (15) 

In Figure [3l we see that the formula (|T5l) is in the complete agreement with the simulated 
results in the case of d = 2. 

The uniform states divide into two classes especially under an external perturbation 
which makes a vehicle accidentally slow down. One is a class of the states recovering 
their uniform configurations, and the other is that of the states never recovering them. 
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Figure 5. Schematic picture of the case that a free vehicle reduces its intention — 1 
approaching a cluster, and finally comes to be in the cluster. The vehicle stops with a 
headway to the cluster ahead, which is estimated from the OV function. 



In this paper, we call such a state unstable against external perturbations a metastatic 
state in accordance with customary practice. Figure H] shows that, under perturbation, 
the flux of a metastable state is decreasing and tends to that of a jam state. (Note 
that it also implies that there is not such a long-lived metastable state as observed in 
[T5].) Consequently, in the region of density where two states coexist, the uniform states 
appear as a metastable, higher-flux branch. 



3.4- The critical point of phase transition 

The flux of traffic flow increases in proportion to the density of vehicles while the 
density is small. However, as the density becomes bigger, close-range interaction between 
vehicles makes a wide, strong correlation over them, and consequently gives rise to a 
jam. Then, there appears a turning point at which the flux declines for the first time. 
Around that point (the so-called critical point), the states of traffic flow bifurcates into 
a stable branch and a metastable branch, and moreover phase transition occurs between 
them. 

In order to estimate the critical point, we consider that the jam line should be 
expressed by 

Q= ^1 (Pmax-P), (16) 

Pmax Pc 

where pmax and pc denote respectively the density at which flux vanishes and that of the 
critical point. From (fT6|) and Q = pv, we have 

l^ = i^^(l-.) + l^.. (17) 

P Pmax Pc 

Equation flTTI) suggests that the spatial pattern divides into two kinds, i.e., clustering 
{v = 0) and free flow (f = 1). Then, since the vehicles move at velocity 1 or in the 
present model, v just indicates the ratio of those in free flow, and moreover total average 
headway (Ax) = (1 — p)/p is calculated from the average headway of the clustered 
vehicles (Axj) = (1 — Pmax)/Pmax and that of the vehicles in free flow {Axp) = (1— pc)/pc. 
These two values reflect a macroscopic property of the SOV model with the OV function 
(fTOj) . but we should estimate them from a microscopic viewpoint. In what follows, we 
consider the case of d = 2, following [22] . 
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Figure 6. The theoretical curve of the maximum density pmax (thick line) at which 
the flux vanishes with the corresponding numerical results (cross). They have perfect 
agreement. We also show a corresponding curve of the original OV model (thin line) 
by use of (^1]) . 



First, we think of pmax as the hmit density at which all free-flow domains of the 
roadway close up and no vehicle can move then. Let us consider the situation illustrated 
in Figure [5] that a free vehicle with its intention 1 is going into a cluster. Then, taking 
the time t = when the headway firstly gets equal to d, the intention decreases as 
f* = (1 — a)*. Since it gives the probability of the vehicle moving at t, the average 
headway, while in cluster, amounts to 



{Axj)=Y[{l-vl) 



t=i 



r7?4(0, 1 - a)%(0, 1 - a)t93(0, 1 - a) 



1/6 



(18) 



(19) 



2(l-a)V4 

where 'dkiu^q) {k = 1,2,3,4) are the elliptic theta functions. Consequently, we obtain 
the density of clustering vehicles 

1 



Pn 



(20) 
2, and 



l + (Axj)' 

as a function of the sensitivity parameter a. Note that < (Axj) < 1 since d 
(Axj) is equivalent to the probability of Aa;* taking the value of 1. 

Figure [6] shows the graph of Pmax, and it has a perfect agreement with the numerical 
results read off from Figure [3l In [22], they give the explicit formula of Axj, the uniform 
headway with which vehicles stop in a jam. It reads 

Ax, = rf-^, (21) 
2a 

where d = 2 and f^ax = 1 in the present case, and a ~ 1.59. By use of (12T1) . we 
illustrate the corresponding graph in Fig. E] as well. Since the original OV model does 
not incorporate a hard-core exclusion rule, the sensitivity parameter a is limited in the 
scope of Ax J > so as to avoid any collision. In contrast, the maximum density of the 
SOV model is retained, due to that rule, not to diverge within < a < 1. 

Next, in order to estimate (Axp) we consider that two vehicles in the front of 
a cluster are getting out of it as illustrated in Figure [71 Then, as described above. 
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Figure 7. Schematic picture of tlie situation that two adjacent vehicles in the front 
of a cluster recover their intention and get out of the cluster. We estimate the headway 
with which the two vehicles finally come to move free. 



there occur two cases since d = 2; = 1 with probabihty (Axj) and = 

with probabihty 1 — (Axj), where we again take the time t = when Ax- becomes 
d. Corresponding to Ax° = and 1, we describe (Axp) respectively as {Axp)^ and 
{Axf)i, and accordingly our main result is expressed as follows: 

where 

{Axf) = {Axf)i{Axj) + (Axi.)o(l - {Axj)). (23) 

Let T denote the interval of time for the front vehicle to get out of the cluster. 
Provided that clusters are large enough to take approximately f° = and that the 
second vehicle leaving the cluster maintains a headway of at least d, we conclude that 
= 0, t;°+i(r) = 1 - (1 - ay, vi = l-{l- ay, and vj+^ir) = 1 - (1 - a)^+* from ([13]). 
Note that r is a stochastic valuable, and hence f and vjj^i are dependent on r. 

In the case of Ax^ = 1: Since vj indicates the probability of moving ahead at one 
site, the probability of r = t amounts to 

Pi{r = t) = viliil ~ vt) (24) 

s=l 

and moreover, the distance which the two make in free flow amounts to 

1 — a 



t) - v] 



a 



vlAr). (25) 



t=i 

Consequently, we have {Axf)i in a convenient form for computation: 

°° 1 - fl 

{AxF)^ =d + Y: —-vl,{r)Pi{r) (26) 

T=l 



_ MO^VTZR (27) 
2(l-a)V8 • 

In the case of Ax^ = 0: The probability of r = t amounts to 

Po{r = t) = vlY.[vt n (28) 

s=l r=l,^s 
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Figure 8. The theoretical curve of the critical density pc (thick line) at which the 
flux bifurcates into a metastable branch and a stable jam branch (i.e. hysteresis), with 
the corresponding numerical results (cross). They also have a perfect agreement. The 
corresponding curve (thin line) of the original OV model is illustrated by use of pop . 



Consequently, from f|T3l) we have 

oo 1 _ ^ 

(Ax^)o = d+Y: [ <i(T)Po(r)] , (29) 



.1- a 



and finally (Axp) is formulated as a function of the sensitivity parameter a. 

The formula of Axp in the corresponding case of the OV model, given in [22], reads 



Ax^ = d+^-^, (30) 

where d = 2, f^ax = 1, and cr = 1.59. Figure [8] shows the critical density pc versus the 
sensitivity parameter a. It also has a perfect agreement with the numerical results read 
out from Fig. [3l By use of fl30|) . we illustrate the corresponding graph in Fig. [8] as well. 
We find that the two theoretical curves present a qualitative agreement, while there are 
some quantitative differences due to the choice of unit car size. 



4. Summary and Conclusion 

The Stochastic Optimal Velocity model was introduced in the preceding paper [15] 
as a stochastic cellular automaton model extending two exactly solvable models (the 
Asymmetric Simple Exclusion Process and the Zero Range Process). Moreover, since 
it has the same formulation as the Optimal Velocity model, the SOV model can be 
regarded stochastic extension of the OV model. 

In the present paper, we take a step function as the OV function in order to 
investigate an elementary property of the SOV model. In previous papers [221 123] . 
the original OV model with this OV function are studied in detail and some analytical 
results are given. Accordingly, we first make a qualitative comparison between them in 
terms of the motion of each vehicle, and then we see a similar hysteresis loop in the 
velocity-headway diagram (phase space) as long as the density is low. That is, as far as 
each vehicle's motion is concerned, the vehicles of the two models show a similar motion 
within a low-density region. 
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However, as the density of vehicles grows large, the structural difference between 
ordinary differential equations and cellular automata, as well as the presence or 
absence of randomness, causes crucial differences especially in the fundamental diagram. 
As discussed in [6], the metastability observed in the OV model is derived from 
the instability of a uniform/homogeneous flow against perturbations. In general, 
randomness introduced in traffic models tends to eliminate metastable states as well 
as unstable ones. Nevertheless, the metastable states appearing in the fundamental 
diagram of the SOV model fend off the disturbance of randomness with a trick, i.e., the 
special choice of the OV function and the configuration. Since that choice allows vehicles 
to move with probability 1, the uniform flow continues to be stable while the uniform 
headway takes the value not less than the discontinuous point of the OV function. That 
simple consideration leads to the formula of the highest-fiux density, i.e., the end-point 
of a metastable branch. However, all the uniform flows are not entirely stable, but 
higher- flux states of them are metastable, i.e., unstable against external perturbations 
as shown in Fig. HI It is one of the most important things to evaluate the critical point 
of density where traffic flow becomes unstable and clustering of vehicles sets in. 

In [7j, they gave the fundamental diagram of the OV model with a practical OV 
function and discussed the stability of uniform flow. Then, they showed the region 
of density where uniform flows become unstable and gave the flux of jammed vehicles 
formulated analytically in the fundamental diagram. In contrast with the original OV 
model, the SOV model incorporates a hard-core exclusion rule, and is thus collision-free 
under any configuration of vehicles. Moreover, as the intention (or average velocity) 
is suppressed since density is high, the vehicle's motion is controlled by randomness as 
well as the hard-core exclusion rule. These aspects give the SOV model a fundamental 
diagram which is obvious different from that of the OV model. 

We point out some apparent differences from the OV model attributed to the effect 
of hard-core exclusion rule and randomness. In the fundamental diagram of the SOV 
model, there is a metastable branch with positive slope (equal to the maximum velocity 
1), and the jam line vanishes at a density less than 1. The simulated results reveal 
that the fundamental diagram of the SOV model with a step OV function consists of 
two line segments with which it has exactly the shape of inversed-lambda. From the 
diagrams, we conclude that the vehicles move approximately with either headway of 
free- flow or that of jammed vehicles, and our issue is thereby reduced to estimation 
of two specific headways. It can be done entirely by probabilistic calculations, and 
consequently the whole fundamental diagram (obviously including the critical point) 
is successfully formulated as a function of the sensitivity parameter. These analytical 
results are attributed to the simplicity of the SOV model and the choice of OV function. 

Further studies on the SOV model under open boundary condition and on the 
multi- velocity version of SOV model will be given in subsequent publications [21]. 
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